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Abstract 

Vesicles under a shear flow exhibit a tank-treading motion of their membrane, 
while their long axis points with an angle < ? with respect to the shear stress if the 
viscosity contrast between the interior and the exterior is not large enough. Above a 
certain viscosity contrast, the vesicle undergoes a tumbling bifurcation, a bifurcation 
which is known for red blood cells. We have recently presented the full numerical 
analysis of this transition. In this paper, we introduce an analytical model that has 
the advantage of being both simple enough and capturing the essential features found 
numerically. The model is based on general considerations and does not resort to the 
explicit computation of the full hydrodynamic field inside and outside the vesicle. 

PACS numbers: 87.16.Dg, 47.60.+i, 87.17.Jj 

1 Introduction 

Vesicles are closed membranes, which are suspended in an aqueous solution. They rep- 
resent an attractive biomimetic system, which has revealed several interesting static and 
dynamical features that bear a strong resemblance with some behavior of real cells. Among 
these features we can cite equilibrium shapes Q] revealing forms similar to red blood cells, 
and tumbling known for these cells |2j. It is known that red cells, like vesicles (Hj, orient 
themselves at a given angle with respect to the shear flow at high hematocrit (high enough 
concentration of red cells), while at low hematocrit (where cells behave as being individ- 
ual) both in vitro and in vivo observations reveal a tumbling motion, where the long axis 
of the red cell rotates in a periodic fashion. It has been recognized for a long time that the 
viscosity ratio between the internal fluid and the ambient one is a decisive factor (the more 
viscous is the internal fluid in comparison to the external one, the easiest is the tumbling). 
Another relevant ingredient is the swelling ratio: a flatten out cell would tumble more 
easily than a swollen one. Several attempts in understanding the tumbling transition have 
been made in the litterature, the most prominent one is the work of Keller and Skalak 
[I]. This work uses the solution of the hydrodynamical equations in the Stokes regime 
(inertial effects are negligibly small for biological blood transport) around an ellipsoid 
which involve quite complex expressions. But still several assumptions had to be made in 
order to solve the problem. Recently, a full numerical analysis has been presented j2] and 
provided the boundaries in the parameter space (basically the viscosity contrast, and the 
swelling ratio) separating the regions of tumbling and those of tank-treading. That work 



focused on vesicles that correspond to a simplified model of red blood cells, and especially 
ignored the elastic properties of the membrane, a fact which though turns out to lead to 
some interesting qualitative changes, will not be accounted for here either. 

Due to the interplay of several effects in the tumbling transition, it is highly desirable 
to have at our disposal an analytical theory, which, on the one hand, should reproduce 
the basic essential features of the tumbling transition, and on the other, should be simple 
enough in order to shed light on the various competing phenomena leading to tumbling. 
It is the main aim of the present paper to deal with this question. 

The present theory bypasses the tedious computation of the velocity field around the 
ellipsoid (Note that in the general case, no analytical solution of the Stokes flow is known), 
and is based on the assumption that the forces acting on each piece of the vesicle membrane 
are proportional to the actual relative velocity at the membrane with respect to the applied 
flow. It follows from our study that simple enough notions account remarkably well for 
many features and render each effect transparent. In addition, this work offers a promising 
basis for more elaborate models, including, for example, the effect of membrane stretching 
or shear elasticity. 

The scheme of this paper is as follows. In section |21 we present the basic ingredients of 
the model. Section |3] is devoted to the derivation of the dynamical equation that governs 
the motion of the vesicle. This part is based on a torque balance. Section |1] presents 
a complementary ingredient that serves to put the evolution equation in a closed form. 
This is based on an energy balance between the inner fluid of the vesicle and the work 
provided by the ambient one. The main outcomes of the analytical theory together with 
their confrontation with the full numerical analysis are presented in section H3 

2 Basic ingredients of the model 

• As in [I], we will assume the shape of the vesicle to be an undeformable ellipse, with 
the long and the short axes denoted by a and b, respectively. It will be recognized 
that the theory can be used for arbitrary prescribed shapes. For definitnessll we 
shall, however, specialize our discussion to an elliptical shape. The enclosed area is 
denoted by S, and the perimeter by P. The fluid embedded into the vesicle has a 
viscosity fii n and the ambient one fi ou t- t is the swelling ratio of the vesicle defined 
in 2D as t = AttS/P 2 . For a circle r = 1 and it is smaller than one otherwise. 

• The vesicle is subjected to a linear shear flow (v x = jy, v y = 0) where 7 is the shear 
rate. Since the vesicles of interest have a fluid membrane, each material point on 
the membrane will be transported by the flow, so that the membrane moves in a 
tank-tread fashion. 

Let us make a remark, which will prove to be useful later in this paper. A simple 
shear flow characterized by the (2x2) shear rate matrix 
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Figure 1: The different frames involved in the model 

can be decomposed into two parts: a symmetric one given by 

7 /2 \ 

7/2 ) 

and an antisymmetric one given by 

7/2 \ 

-7/2 ) 

As shown on Fig. |5J the antisymmetric part provides a rigid-like clockwise rotation of 
the vesicle (R), while the symmetric part corresponds to an elongational (or strain) 
flow, which tends to orient the vesicle along ? (E). 

Our calculation is based on the following two properties of the Stokes equations: 

• Due to the linearity of the Stokes equations, the superposition principle for given 
boundary conditions applies: the velocity field around a vesicle subjected to a tank- 
treading and a tumbling motion in a simple shear flow is the sum of the velocity 
fields obtained for the three following configurations (see Fig. |3J): 

— A simple shear flow acting on a rigid body fixed in the flow at a constant 
orientation angle ip with a fluid velocity equal to zero on the contour of the 
vesicle. 
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Figure 2: Decomposition of the shear flow (S) in a rotational part (R) and an elongational 
part (E). 

— The flow created by a rigid elliptic body rotating at a rotation velocity ^ in a 
quiescent fluid. 

— The flow created by an elliptic body subjected to a tank-treading motion of its 
contour and fixed at a constant orientation t/j in a quiescent fluid. 

• The second ingredient, which follows from the previous one, is an extension of a 
general result valid in Stokes flows for a solid which is in relative motion at a ve- 
locity V with respect to the surrounding fluid. The drag force on the solid scales 
as Fdrag = A* A V, where A is (a drag coefficient) function of the geometry of the 
body. 1 There is a linear relation between the force and the relative velocity of 
the body with respect to the applied flow. We view the elliptic contour as being 
represented by adjacent segments. The key hypothesis of the following analysis is to 
apply this property, ie. the linearity between forces and relative velocities, on each 
segment of the membrane. Let us make some important comments about the mean- 
ing of this assumption. The external force applied on an elementary segment of the 
membrane is provided, on the one hand, by the flow imposed externally, and, on the 
other hand, by the back-flow due to the presence of the vesicle. This retroaction of 
the vesicle on the applied flow is a complex piece of the study and an exact determi- 
nation of its effect requires sophisticated numerical treatments such as the Boundary 
Integral method ([El El)- In our model, the basic assumption stated above takes into 
account this complex interaction in an effective manner: the effect of the back-flow 
is included in the coefficient A, which links the effective force to the relative velocity 
of the segment with respect to the applied flow. In the framework of our model, 

1 More precisely, the Stokes force exerted on a solid of typical length L in a translational motion at 
speed U, in a quiescent newtonian fluid of viscosity scales indeed as 

F « fiUL 

More formally, we can write a linear relation between the force and the velocity: 

F t = -nAtjUj (1) 

Aij is a tensor which is symmetrical for a newtonian fluid, and in a specific frame linked to the solid, one 
can write : 

Fi = -fiXiUi (2) 
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Figure 3: Decomposition of the velocity field around the vesicle subjected to a tumbling 
motion and a tank-treading motion of its membrane in a simple shear flow. 

this coefficient is chosen to be independant on the particular elementary segment 
considered. This is reminiscent of a "mean-field" like approximation. This coeffi- 
cient is also a priori non isotropic, i.e. takes different values depending on wether we 
consider the normal direction or the transverse direction of the elementary segment 
considered. These two values will be denoted as A* and A n and their determination 
will be discussed in section To some extent this model is akin to the Rouse model 
for polymer rods jSj, where hydrodynamical interactions between adjacent pieces are 
ignored. 



3 Mechanical equilibrium for the vesicle in the shear flow 

As stated above, the main idea is to use the linear generalized Stokes law at the local 
level of each segment of the contour, and to compute the torque associated with the force. 
Since we shall decompose the velocity field into an applied shear, a tank-treading motion, 
and a tumbling one, we shall have to deal with three types of forces separately. Once each 
torque is evaluated, we sum up the three contributions, and set the resultant to be zero, 
owing to absence of inertia. Once the expressions of the forces are specified, the remaining 
pieces of the work are purely algebraic with some specific integrals involving the geometry 
of the vesicle. 

Since the Stokes law relating the force to the relative velocity is local, we find it 
convenient to first, write it in the frame linked to the vesicle, and then, to express the 
torque elements in the laboratory frame for ease of computations. We refer to Fig. ^for 
the different frames used here. The laboratory frame has the basis denoted as (e x ,e y ). 
The rotating frame, which is linked to the principal axes of the ellipse, is specified by 
its basis denoted (e' x ,e' y ). The local frame associated to an elementary segment on the 
elliptic contour is specified by (e t ,e n ). 

In the local frame (et,e n ), the components of the local force (or drag) applied on the 
segment per unit length in the transverse direction can be expressed as functions of the 
relative velocities (Vt,V n ) exerted upon a membrane element: 
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dF t = -n out \ t V t dl ^ 
dF n = — /j, ou t X n V n dl 

where At and X n are phenomenological parameters of the model associated respectively to 
the transverse and the normal motions to the segment dl. Xt and X n have positive values 
with the choice of Eqs. 01 They have a dimension of the inverse of a length. The crux of 
the analysis is to decompose the local velocity in three pieces as stated above, and evaluate 
various torques. 

• (1): Torque of the force acting on the body in a simple shear flow 

In the laboratory frame (e x ,e y ), the velocity field of a simple shear flow takes the 
form: 



v x = jy 

Vy=0 

7 is the shear rate, which fixes the time-scale of the flow. Written in the local 
coordinate system (et,e n ), the relative velocity reads: 



Vshean = ~1 U^-^t 
Vshear n = 7 2/6 x .e n 



Using (J2J , we determine the associated forces denoted as dF s y ieart and dF s h earn , from 
which the torque is computed as: 



M sh ear = j> * X dF shear (4) 

Using the coordinates linked with the natural axes of the ellipse (for a convenient 
calculation), we easily find: 



M s hear ~ HoutJ 



L 2 s i + L 2 s2 L 2 sl - I? 



+ 



,s2 



cos(2?/>) 



(5) 



with the convention that a positive torque corresponds to a clockwise rotation (see 
Fig. L S 2 and L s i are elliptic integrals: 



L 2 s2 = [\ t f(a, b) - X n f(b, a)] + (£) 2 [(X t g(a, b) + \ n f(b, a)] 
L 2 sl = [X t f(a, b) - X n f(b, a)] + (f ) 2 [(X t g(b, a) + X n f(b, a)] 

with 

f(a, b) = § c ( . )V ^^ ai ,, di 

g (a, b) = § c (a )2y ,? + (b )2x ,2 dl 



According to the linear decomposition of Fig. [21 we can identify from (jSJ the torque 
associated with the rotational part of the flow: 
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M rot = Houtl g ^> 

and the torque associated with the elongational part of the flow: 

L 2 - L 2 

M e i ong = /W7 sl 2 s2 cos(2V>) (7) 
(2) Torque of the force acting on a rigid ellipse with a rotation speed ^ 

The tumbling velocity at a position r of the membrane is given by: 

Vtumble = w x r (8) 

where w = ^e z is the instantaneous angular velocity of the vesicle 
and its components in the local frame are 



Vtumblet dt^y' ^ dt " x ') 

dt *V » Gx 



V tumblen = (+*'f e y , - y'f e 



We use © to determine the force and then we compute the torque as: 

Mt um ble = j c rX dF tumble- 



IC 

This yields, after elementary integration, to: 

^^/^ / 2 2 \ / \ 

Mtumble = ^out~l^{L sl + L s2 ) (9) 

(3) Torque of the force acting on the ellipsoid related to the tank-treading motion 
The tank-treading velocity is tangential to the membrane: 

Vtank = Vtank^t (10) 

and the force is simply: 

Ftank = —fJ-out^tVtank^t (H) 

The associated torque is: 



M tank = j> r x dF tank . 

The integration provides us with: 

Mtank = -HoutVtank(Li + L 2 ) (12) 
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where 



In = A t / c (V 2 )/y[(f)V2 + (f)V 2 
L2 = Xt$ c (ly' 2 )/J\(V 2 x' 2 + &y' 2 



dl 
dl 



Because the inertial effects are small (and thus neglected), the sum of the three torques 
must vanish. Summing up the three contributions (Eqs. ©, © an d (|12j0 . one finds the 
evolution equation for the angular velocity of the vesicle: 

d$> 

— — = LO rot + Ll> c + belong, (13) 

where we have defined the three quantities on the r.h.s. of (|13(l as 

1 

Urot = — ^ ( 14 ) 

where oj ro t represents the rotational velocity (i.e. a torque in the Stokes framework), 
arising from the rotational part of the flow and is responsible for the global rotation of 
the shape. 

W c = Vtank-pi 1 — TT (15) 

L sl + L s2 

lo c is the contribution of the tank-treading motion of the membrane to the effective angular 
velocity ^ of the vesicle. 



2 I 2 sl + L 2 s2 



belong = , j€ COS ( 2 ^) ( 16 ) 



belong is the effective elongational velocity which represents the main contribution of the 
elongational flow and tends to orient the vesicle along a direction making an angle if) = ? 
with respect to the applied flow. 

It is interesting to note at this point that for a sphere, the various integrals can easily 
be computed 

Li + L 2 = 2ir\ t a 2 ; I 2 sl + L 2 s2 = 2A t vra 3 

Reporting into ()14|) - (|15|) . and requiring in Q13|) that d^//dt = (since for a sphere the 
contribution to tumbling vanishes 2 ), we obtain: 

Vtank = T^a (17) 

2 Distinguishing between tumbling and tank-treading for a sphere might seem a bit confusing. The case 
of a sphere is degenerate, since one can view the dynamics as being of pure tank-treading or pure tumbling 
nature. Indeed requiring either that = 0, as we did here, or Vtank = , provides the same velocity 
along the contour. For continuity reasons with the case where there is a slight deviation from a sphere, we 
interpret the motion under question as being of tank-treading type. 
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This is the expected value of the tank-treading velocity in the case of a sphere with a 
radius a. Interestingly, this result holds whatever the prescription for the parameters A n 
and \t- 

Inspection of eqs (|13|) - (|16|) reveals, in particular, that for a rigid membrane where 
Vtank = 0, no stationnary solution is possible: a rigid elliptic body should always tumble, 
as one expects. This can be interpreted by the fact that the rotational velocity | iv ro t | is 
always bigger than the elongational velocity | uJ e i on g !• If allowance is made for a tank- 
treading motion (due to the membrane fluidity and the finite viscosity of the internal 
liquid) then uj c ^ 0. Equations (fT3)h ((11]) and (fT3|) show indeed that the tank-treading 
motion described by the velocity co c results in an effective reduction of the global rotation 
^rot, provided that the tank-treading velocity Vtank has a positive value. A stationary 
(non tumbling) tank-treading motion of the shape is thus possible if the velocity | u> e iong |, 
representing the elongational part of the flow, can balance the effective tumbling velocity 
| Lu ro t + uj c |. This can occur for sufficiently high values of the tank-treading velocity. 
Hitherto, the tank-treading velocity has been introduced as a phenomenological quantity, 
and it must be computed independently. This step is necessary in order to have an 
evolution equation in a closed form. The tank-treading velocity is clearly limited by the 
viscous friction of the internal fluid, and this piece of information must be evoked in order 
to complete the analysis. 



4 Determination of the tank-treading velocity 

Following [3], the tank-treading velocity Vtank can be determined by considering the energy 
dissipated in the system. The energy injected by the flow is dissipated by viscous friction 
in the fluid inside the vesicle. The energy rate (or power) provided by the fluid to an 
elementary segment in the laboratory frame is equal to dF.Vt. 

The velocity Vt at a point r of the membrane can be written in the same frame as: 

V t = V tank e t + w x r (18) 

where w = ^ez is the instantaneous angular velocity of the vesicle. 

An elementary force dF acting on an element dl of the membrane can be decomposed 
according to the previous section as: 

dF = dF s hear + dF tum ble + dF tan k 

The sum of the torques applied on the vesicle is equal to zero, entailing: 

j> dF.(w x r)dl = <j> w.(r x dF)dl = (19) 

Hence, only the first contribution of the velocity in (|18|) matters. The total power provided 
by the flow has the following contributions: 

Etot = E s h ear + Etumble + E t ank (20) 

• The contribution from the simple shear flow is E s h ear = § dF s h e ar-^4anfc e t 
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and upon integration on the contour, we find: 

L\ + L2 L\ — L2 



E. 



shear 



V-outVtankl 



+ 



cos(2^) 



(21) 



where the lengths L\ and L2 have been defined previously (easll2[). Following the 
spirit of the last section, we write E s h ear = E rot + E e i ong in order to identify the 
contributions from the rotational part of the shear flow 



E, 



rot 



V-outVtankl 



Li + L2 



and the elongational part of the shear flow 

Eeiong = ^outVtankl^—^^ COs(2V>) 

This decomposition will be useful in the discussion of the results. 

The contribution from the tumbling motion is E tum bi e = § dFtumbie-V^^et, yield- 
ing: 



dip 

Etumble = V-outVtank-^\L\ + L2) 

• The contribution from the tank-treading motion is Et an k = § 
upon integration one finds: 

Etank = — ^out^tankE' 

where P' = § c X t dl = \ t P 
By using the above results, the total power (|2()|) takes the form: 

Etot — H-out 

{oiV tank + (3Vtank) 



(22) 

dF taak .V tank e t , and 

(23) 



where 



a 



(L 1 + L 2 ) 2 
L 2 S \ + L 2 S 2 



P' 



(24) 



(25) 



/3 = 7 



Lo — Li 



L\ + L 2 L 2 s2 



2 ^^(2*) 



2 L s2 + L sl 



(26) 



In the particular case of a circular shape, the total power provided by the external fluid 
to the internal one can easily be determined: indeed, we have (L\ + L2) 2 = 4\ 2 (ttci 2 ) 2 ; 

+ L 2 2 = 2Aj7ra 3 ; P' = 2Af7ra, and L\ = L2, L s \ = L S 2- This implies that both 
coefficients a and (3 vanish, and so does the total power. This result is conforting since 
inside a sphere the fluid executes a rigid-like rotation (there is no dissipation) and thus no 
energy can be transferred. It is only when the shape deviates from a circle (or a sphere in 
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3D) that dissipation is permissible. Note that we arrived at this result before using any 
information about dissipation in the enclosed fluid, and this points to a consistency of the 
model. 

The energy dissipated by viscous friction in the volume of the vesicle is of the form: 

1 f.dVt dV U2 , 

In general, we have to determine the velocity field, which satifies the Stokes equations 
inside the vesicle and subjected to boundary conditions at the surface of the ellipse. Our 
aim is not to determine the velocity field exactly, which is not an easy task in general (and 
an exact result is the exception rather than the rule) . Rather we wish to capture the main 
ingredients and remain within a heuristic analysis. For that purpose, it will be sufficient 
to make use of an approximate solution inferred from simple considerations based on the 
result relative to a spherical shape. In order to anticipate the main ingredient, we shall 
take the case of a slightly deformed circle as a reference in order to serve as a guide for our 
reasoning. Consider e = to be small. The following velocity field fulfills the prescribed 
conditions (i.e. to be a solution of the Stokes equations in the inner domain of the vesicle): 

Vx' = Vtank\ 
Vyi = —Vt an kX 

It must be noted that despite the fact that the velocity is not exactly constant along 
the contour, the velocity remains colinear to the tangent at the membrane 3 and this 
continues to represent a reasonable approximation. Let us estimate the energy dissipated 
in the vesicle. This is given by : 

e = fi in a'V t l nk , (28) 

where a' is a constant depending on a and b: a' = ^irab(\ — ^) 2 in the present case. 
Using (|2i|) and (|28|) . we arrive at: 

Vtank = — f f l„ , 7cos(2*) (29) 

where 



.1,-16. 
a = — 7rao(— - -=) (30) 



2 K b 



{Li+L 2 f pl 
L 2 S \ + L 2 s2 



(31) 



3 Other prescriptions for the flow could have been used. In particular: 



V y l - ~Vtank^ 



This flow ensures a constant value for the tank-treading velocity along the contour but the velocity is not 
colinear to the tangential direction of the contour. 
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L 2 - Lg Li + L 2 L 2 s2 - L 2 sl 

h ~ ~ (32) 

As could be anticipated, the tank-treading velocity is directly proportional to 7 which 
fixes the time-scale of the imposed flow. We also note that the tank-treading velocity 
enjoys the same symmetry as the elongational flow does: it vanishes for i[) = ±? and is 
maximal for ip = 0. 

Some remarks are in order. As explained in section 13 a shear flow can always be 
split into a rotational part and an elongational one (see Fig. |2J). In order to understand 
the origin of the tank-treading motion, it is appropriate to specify the role of both the 
elongational and the rotational components of the flow. 

On the one hand, for a purely rotational flow, eq. lf24*j) shows that E to t = ^outotV 2 ank 
since (3 = 0. Equating and (|2*S|) leads to the condition Vtank = 0, provided that 
the shape is not circular. This corresponds to a global solid-like rotation. On the other 
hand, for a purely elongational flow, eqs. (|2*1|) and (f2*%|) lead to Vtank ~ cos (2^;). A 
non-zero tank-treading velocity is possible with the proviso that the orientation angle is 
different from ip = ? (cos(2^) 7^ 0). The torque applied on the vesicle arising from the 
elongational flow is (see eqJTJ) M e i ong ~ cos(2^). An inspection of the balance of the 
torques for the elongational flow, as done in section |3J leads to the dynamical equation 
£jt ~ cos(2ip). ip = j is thus, the only steady equilibrium position, with a tank-treading 
velocity equal to zero. Hence, the existence of a tank-treading motion of the membrane is 
only a consequence of the coupling between the rotational and the elongational part of the 
flow. The total effect of the shear can be interpreted as follows: the rotational part tends 
to push the orientation angle of the vesicle axis towards lower values than tp = ? . As soon 
as this is achieved the vesicle acquires a non-zero tank-treading velocity since there, the 
elongational part enters into action (see eq. (|16|0 4 

5 Dynamical equation for the orientation angle 

Plugging eq. (|29|) into eq. we can express explicitly uj c (which involves the tank- 

treading motion of the vesicle) and this leads to the general dynamical equation for the 
orientation angle ip of the vesicle: 

= A + B cos(2iP) (33) 

dt 

where 

A = -l (34) 

4 Note that the tank-treading velocity is the result of an energy balance which involves the coupling 
between the rotational and the elongational parts of the flow. Since energetic quantities are quadratic func- 
tions of the velocity field, the tank-treading velocity is not simply the sum of the tank-treading motions 
associated respectively to the elongational component and the rotational component considered indepen- 
dently. Such a summation would result in a vanishing tank-treading velocity 
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D 



1 
2 



L 2 



s2 



L 2 



.si 



+ (L 2 +£l) 



(L 2 - L x ) - (Li + L 2 



.«2 



(L l + L 2 y-(P' + -^-a')(L^ + L% 



(35) 



A purely tank-treading motion corresponds to the situation where the inclination angle 
is constant. This is expressed by = which implies the condition: — -^ < 1. This 
constraint leads, in particular, to a condition on the viscosity ratio between the inner and 
the outer fluid: 



fJ'in 1 



L 2 



(L 1 + L 2 )(l+ r2 



si 



T2 I j2 
^s2 L s2 



P' 



(36) 



This is the general condition which can be tabulated numerically, provided that the two 
drag coefficients At and A n are known, which is exactly the case for several shapes (disks, 
ellipses, spheres). Thus, the condition relates uniquely the viscosity contrast to geometrical 
quantities which are functions of the swelling ratio. In order to gain more insight towards 
an analytical progress we can explore the situation of a small deformation around the 
spherical shape. This proves to be sufficient to capture the essential features. For that 
purpose, we set e = b — a and treat e as a small parameter. The integrals L\,L%... that 
enter in ()36|) can be evaluated explicitly, so that the critical condition for tumbling is 
expressed in a simple form in terms of the viscosity ratio and the swelling ratio r: 



/'if 



5Xta 



Vout 1 - T 

and the dissipation rate scales as: 



(37) 



(38) 



This law for the dissipation rate is in a good agreement with previous numerical results 
(see jS]). We find here that the closer is the shape to a sphere, the more difficult does 
tumbling occur; the viscosity contrast for tumbling diverges as 1/(1 — r). 



6 Quantitative and qualitative analyses of the model 

We have seen that the model presented here captures the essential features and sheds 
light on the various competing effects that fix the tank-treading and tumbling motions. 
We may ask the question whether the model can be made more quantitative. As stated 
before, the model requires the introduction of two drag parameters Aj and A n which 
are the proportionality constants relating the force and velocity in the normal and the 
tangential directions. Let us recall that these two drag coefficients describe the effect of 
the hydrodynamic interaction felt by a membrane element. This involves the geometry 
of the shape around a given element of membrane and these drag parameters are a priori 
function of the position of the element considered on the contour. If one wishes to go 
beyond a qualitative discussion we must determine these two drag parameters, which can 
be made in general only numerically. For a sphere with a radius a moving in a Stokes 
flow, it is known that on a local segment of the spherical contour we have: At = A„, = 3/2a 
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(see [0]). Instead of evaluating the exact values of these parameters for each elementary 
elements of the membrane, we shall rather estimate them from the best fit with the full 
numerical simulations obtained previously(|Hj). Note that each membrane portion can be 
approximated locally as an arc of a circle. Thus, as in the case of a sphere we chose equal 
values for the drag coefficients in the normal and in the transverse directions with respect 
to the contour: Aj = A n . We consider now a vesicle with an aspect ratio r = 0.8 and 
determine several quantities. 

The results have been confronted to the full numerical computation and the values of 
At and A n have been guessed. For various tests made so far, we found that Xt = X n — 4 
per unit length in the transverse direction provide the most reasonable fit. We present on 
Fig. 01the evolution of the equilibrium angle as a function of the viscosity ratio r, which 
corresponds to the stable branch of the saddle-node bifurcation JOJ • The prediction of the 
model qualitatively reproduces the bifurcation branch (this is always the case regardless of 
the chosen parameter), and is fairly in reasonable agreement on the quantitative level. The 
point at which the angle is zero corresponds to the threshold of the tumbling bifurcation. 
This threshold depends on the swelling ratio. Consequently, the two parameters controlling 
the bifurcation are the viscosity contrast and the swelling ratio. Fig. |S] represents the 
boundary between the region of the phase diagram where pure tank-treading motion takes 
place (low r and high r) and that, where the motion is of tumbling type (which is favored 
at large r and small r). The results are compared with the full numerical calculation. It 
is also worthwhile to represent some other physical quantities. Of particular interest are 
the global rotation velocity and the tank-treading velocity (Figs. 16171 and I8I9|) . 

Before concluding some additional comments are worth to mention. Following the 
considerations in section 21 the tank-treading motion is a result of the competition be- 
tween the rotational part of the flow and the elongational component. More precisely, 
the rotational component pushes the vesicle axis away from the elongational main axis 
(t/j < it/ 4), allowing the membrane to acquire a non-zero tank-treading motion. One may 
say that part of the rotation torque is transferred to the tank-treading one. Increasing the 
viscosity of the inner fluid results in a global reduction of the tank-treading velocity since 
the internal dissipation penalizes velocity gradients inside the vesicle (see Fig. EJ). From 
Fig. IHlthe effective tumbling velocity should thus increase, reducing further the value of ij). 
However, thanks to the cos(2?/>) variation of the tank-treading velocity, a new equilibrium 
position can be found at a value of ifi which is a decreasing function of the viscosity ratio. 
In the extreme limit where ip eq = 0, the elongational velocity uj e i ong reaches its maximum 
and can not overcome the rotational velocity on further increase of the internal viscosity: 
the steady-state solution doesn't exist anymore whereby a new dynamical solution takes 
place in the form of tumbling. Fig. 0] illustrates the evolution of the equilibrium angle as 
a function of the viscosity ratio. 

The above discussion was made on the assumption that the swelling ratio is constant. 
The swelling ratio is a measure of the deviation from a spherical shape. For the particular 
case of a circular shape (corresponding to a swelling ratio r = 1), the total velocity arising 
from the rotational part of the flow is completely transferred in the tank-treading motion 
of the membrane so that the effective tumbling velocity | uo ro t + 00 c \— (in reality, as 
commented above, this situation is degenerate and there is no distinction between rigid 
rotation and tank-treading). For a circular shape, the tank-treading velocity is maximal 
and equal to V tan k = 2 a wnere a 1S the radius of the circle. This result has already been 



14 



obtained directly in section 03 

Figures H3 and [7| represent the evolution of the tank-treading velocity and the effective 
tumbling velocity | u ro t + uj c | as a function of the swelling ratio r. As the shape deviates 
from a circular one, the effective velocity responsible for tumbling | uj ro t + tv c | increases 
(Fig. EJ) since the tank-treading velocity decreases (Fig. EJ). This explains that the 
transition to a tumbling regime can be achieved for lower values of the viscosity ratio as 
the swelling ratio decreases. This is indeed what is observed on Fig. |SJ 

7 Conclusion 

We have presented a simple model bypassing the calculation of the Stokes flow. We have 
captured the essential features of the transition tank-treading/tumbling, and have a trans- 
parent view of the various competing phenomena. This work has added a piece to our 
understanding of tumbling. There are several important effects that have been disre- 
garded, however. We have restricted most of our discussion to 2D shapes. In view of the 
result of Keller and Skalakjl] we do not expect a qualitative change when 3D shapes are 
considered, provided that the shape is prescribed. If the shape is free to evolve, the shear 
may induce a shape transformation, like a prolate/oblate transition, and this constitutes 
an important task for future investigations. In addition, we did not include the fact that 
the two monolayers forming the vesicle membrane may slide with respect to each other 
[TT] . In that case one has to include two tank-treading velocities, one for each layer, and 
evaluate the membrane internal dissipation. It will be an interesting point to clarify the 
influence of this fact on tumbling. For biological cells, like red cells, further refinement of 
the model is clearly necessary. For example, red blood cells tumble in vivo in the same 
manner as vesicles do. There is however a notable difference between vesicles and red cells. 
The transition to tumbling depends on the shear rate 2 . This dependence is completely 
absent for vesicles since there is only one time scale I/7 which is imposed by the flow. 
This points to the fact that there should exist a relevant intrinsic time scale for red cells. 
A natural candidate is the elastic (or even viscoelastic) response of the cytoskeleton. A 
natural time scale is [ijG where G is the (2D) shear modulus of the spectrin network 
forming the red cell cytoskeleton, and \x is the membrane viscosity. Available data on G 
[12] and \x ^2] provide us with fi/G of the order of 1CT 2 — I CP 1 seconds, which is not far 
from I/7 in ordinary experiments We are presently using a simple model for elasticity 
in order to analyse the qualitative features of this effect Jl] . 

This work has benefitted from a final support from ONES (Centre National d'Etudes 
Spatiales). 
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Figure 4: Equilibrium angle as a function of the viscosity ratio: the saddle-node bifurcation 
(r = 0.84: model; r = 0.8: simulations) 
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Figure 5: Evolution of the critical viscosity ratio rc as a function of the swelling ratio r 
for a = 1. 
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Figure 6: Evolution of the effective tumbling velocity (| uj ro t + lo c | /7a) as a function of 
the swelling ratio r, r = = 2 




Figure 7: Evolution of the tank-treading velocity Vt an k /7a as a function of the swelling 
ratio r 
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Figure 8: Evolution of the effective tumbling velocity ((| uj ro t + uj c \)h a ) as a function of 
the viscosity ratio r, the swelling ratio is equal to r = 0.84 
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Figure 9: Evolution of the tank-treading velocity Vtanklio, as a function of the viscosity 
ratio r, r = 0.84 
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